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ABSTRACT 

Recent data from the Tevatron have revealed that the production rate of prompt 
charmonium is more than an order of magnitude larger than the best theoretical 
predictions of a few years ago. This surprising result can be understood by 
taking into account new production mechanisms that include fragmentation and 
the formation of charmonium from color-octet cc pairs. 



1. Introduction 



The typical high energy physics conference these days includes talk after talk 
showing remarkable agreement between experiment and theory. There is an occasional 
two-sigma discrepancy, but most such problems will go away if you have the patience 
to wait for better data. However there is one problem where experimental results 
have differed from theoretical predictions by orders of magnitude. This problem is 
the production of charmonium at large transverse momentum at the Tevatron. 



2. Color-singlet Model 

Until recently, the conventional wisdom on the production of heavy quarkonium 
was based primarily on the color-singlet model^ In this model, the cross section 
for producing a charmonium state is proportional to the perturbative cross section 
for producing a color-singlet cc pair with vanishing relative momentum and with 
appropriate angular-momentum quantum numbers: ^Sq for r^^, ^5*1 for J/ip and ip' , 
^Pj for Xcj, etc. The color-singlet model has great predictive power. The cross section 
for producing a quarkonium state in any high energy process is predicted in terms 
of a single nonperturbative parameter for each orbital-angular- momentum multiplet. 
The nonperturbative factor is |i?(0)p for S-wave states, |i?'(0)p for P-wave states, 
etc., where R{r) is the radial wavefunction. For example, the inclusive differential 
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cross sections for producing J/ip and Xcj in the color-singlet model have the form 



+ X) 
dcriXcj + X) 



da{cc{i,^Si)+x) \R4o)\^, 

da{cc{h'Pj)+X) \R'^S^)\'. 



(1) 
(2) 



As the name suggests, the color-singlet model is not a complete theory of quarko- 
nium production derived from QCD. The model ignores relativistic corrections, which 
take into account the nonzero relative velocity v of the quark and antiquark. These 
corrections may be numerically significant, since the average value of v'^ is only about 
1/3 for charmonium and 1/10 for bottomonium. The color-singlet model also assumes 
that a cc pair produced in a color-octet state will never bind to form charmonium. 
This assumption must break down at some level, since a color-octet cc pair can make 
a nonperturbative transition to a color-singlet state by radiating a soft gluon. The 
clearest evidence that the color-singlet model is incomplete comes from radiative 
corrections. In the case of S-waves, these can be calculated consistently within the 
color-singlet model. However, in the case of P-waves, the radiative corrections con- 
tain infrared divergences that cannot be factored into |-R'(0]P- This problem was 
first noted in 1976 in connection with the decays of Xc states,0 but it was solved only 
recently.i The divergence arises from the radiation of a soft gluon from either the 
quark or the antiquark that form the color-singlet ^Pj bound state. The infrared 
divergence can be factored into a matrix element (Og^l^Si)) that is proportional 
to the probability for a pointlike cc pair in a color-octet ^5*1 state to form Xc plus 
anything. Thus, perturbative consistency demands that the formula of the color- 
singlet model be modified to take into account the nonzero probability for a cc pair 
produced in a color-octet state to bind to form charmonium: 



The color-singlet model can be used to predict the production rate of charmonium 
at large transverse momentum in hadron colliders. The first thorough treatment of 
this problem was given by Baier and Riickl in 1981, and their analysis remained the 
conventional wisdom for the next decade.! A ip with large pt can be produced either 
directly, or from a Xcj with large pt that decays via Xcj — > ^7, or by the decay of 
a B hadron with large px- Baier and Riickl assumed that the direct production of 
charmonium is dominated by the parton processes that are lowest order in the QCD 
coupling constant a^. The relevant parton processes that produce cc pairs at large 
Pt are gg ^ cc + g, gq ^ cc + q, gq ^ cc + q, and gg — >• cc -t- g. The cross sections 
da for these processes are all of order a^, but they have different dependences on 
Pt- The only parton process that produces direct ip is gg ^ cc + g, and it gives a 
cross section that has the behavior dd/dp^ ~ l/p|, at large pt- The dominant parton 
process for direct XcJ 99 cc + g, and it gives da /dp^ ~ 1/Pt- Both of these 
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cross sections fall more rapidly with pt than typical jet production cross sections, 
which behave like da/dp^ ~ 1/Pt- The cross section for b quark production has this 
scaling behavior when pt ^ mb. Thus the conventional wisdom was that ip^s at large 
Pt should come predominantly from b quarks, with direct Xcj's being the next most 
important source, and direct ip^s being negligible. This conventional wisdom has been 
completely overthrown by recent experimental data from the Tevatron. 

3. Prompt Charmonium at the Tevatron 

In the 1988-89 run of the Tevatron, the cross section for ip production at large 
Pt was measured by the CDF collaboration. They also measured the cross section 
for Xc at large pr- Assuming the conventional wisdom that V^'s at large pr come 
predominantly from b quarks and from direct Xc'Sj^ they inferred that the fraction 
fb of ip^s that come from b quarks was about 60%.l3 This number was then used to 
determine the 6-quark cross section. Unfortunately, the conventional wisdom proved 
to be wrong. The fraction is actually closer to 15%, and the b quark cross section 
is significantly smaller than the result obtained in the 1988-89 run. 

The breakthrough came in the 1992-93 run of the Tevatron, with the installation of 
a silicon vertex detector at CDF. This device can be used to measure the separation 
between the collision point of the p and p and the point where the ip decays into 
leptons. If a ^ with large px is produced by QCD mechanisms, then the leptons from 
its decay will trace back to the pp collision point and the ip is called prompt. Similarly, 
if a Xc is produced by QCD mechanisms and decays radiatively, the resulting ip is 
also prompt. On the other hand, the ip^s coming from b quarks are not prompt. A 
B hadron with large px will travel a distance on the order of a millimeter before it 
decays weakly. Thus the leptons from the decay of the ip will trace back to a vertex 
with a measureable displacement from the collision point. 

By using the vertex detector to remove the background from b quarks, CDF was 
able to study the QCD production mechanisms directly. They discovered that the 
cross section for prompt ip production at large pt is about 30 times larger than 
predicted.il The prompt cross section for Xci + Xc2 is also larger than predicted by a 
similar factor. The cross section for prompt ip' at large pt is three orders of mag- 
nitude larger than the theoretical prediction. The conventional wisdom on prompt 
charmonium production fails completely when confronted with this data. The data 
can only be explained by new production mechanisms. 

4. Fragmentation 

Braaten and Yuan pointed out in 1993 that the dominant production mechanism 
for charmonium at sufficiently large pt must be fragmentation, the formation of char- 
monium within the jet initiated by a parton with large transverse momentum.! This 



production mechanism had not been included in any previous theoretical predictions 
for pp collisions. The importance of fragmentation should have been realized long 
ago. There are factorization theorems of perturbative QCD dating back to 1980 that 
guarantee that the inclusive production of a hadron at large pt is dominated by 
fragmentation.0 According to these theorems, the asymptotic form of the inclusive 
differential cross section for producing a hadron H with momentum P is 

dcr{H{P) + X) = dz dd{i{P/z) + X) Di^niz) , (4) 

. Jo 

where da is the differential cross section for producing a parton of type i with mo- 
mentum P/z and Di^niz) is a fragmentation function. It gives the probability that 
a jet initiated by parton i includes a hadron H carrying a fraction z of the parton 
momentum. The parton cross sections da can be calculated using perturbation the- 
ory (up to parton distributions for the colliding particles if they are hadrons). All the 
nonperturbative dynamics involved in the formation of the hadron H is contained in 
the fragmentation functions. 

The factorization theorem (^) applies equally well to heavy quarkonium as to 
light hadrons. It is easy to see that the conventional calculations of ip production 
at large px did not include fragmentation contributions. At large pt-, the parton 
cross sections da in have the scaling behavior da/dp\ ~ up to logarithmic 
corrections. However the leading-order cross sections for gg ip + g and gg — Xc.j + 9 
in the color-singlet model behave asymptotically like l/pf^ and l/pf^, respectively. So 
where are the fragmentation contributions? The answer is that they do appear in the 
color-singlet model, but only at higher order in perturbation theory. They appear 
at next-to- leading order for Xc and at next-to-next-to- leading order for ip. These 
higher-order corrections would be prohibitively difficult to calculate in their entirety, 
but fortunately the fragmentation functions can be calculated relatively easily. The 
fragmentation functions that are needed for prompt ip production in pp collisions, 
including g ^ ip, c ^ _th, g ^ Xcj, c — > XcJ, and 7 — have all been calculated 
to leading order in a^.B The fragmentation functions for tp were calculated in the 
color-singlet model. The fragmentation functions for Xcj were calculated using the 
factorization formula (^, which includes a color-octet contribution. 

To calculate the fragmentation contribution to the cross section for prompt ip 
production at large pt in pp collisions, the fragmentation functions Dg^uiz) and 
Dc^h{z) must be folded with parton scattering cross sections for ij ^ g + k and for 
ij — > c+k, c+k, respectively. They must then be convoluted with parton distributions 
fi/p{xi) and fj/p(x2) for the colliding p and p. The largest contribution comes from 
the color-octet term in the fragmentation function for g Xc- The fragmentation 
contribution increases the theoretical prediction by an order of magnitude, bringing 
it to within a factor of 3 of the CDF data.i Given the many sources of uncertainty in 



the calculation, this can be considered reasonable agreement. 

5. Color-octet Mechanism 

While fragmentation may explain the CDF data on prompt ip production, it is 
not enough to explain the data on prompt ip'^s. In this case, the conventional wisdom 
gave predictions that were 3 orders of magnitude smaller than the data. When frag- 
mentation is included, the prediction increases by more than an order of magnitude 
at large pt-, but it remains more than an order of magnitude below the data. The big 
difference between tp and ip' is that the ip signal is fed by the decay of Xc's, but there 
are no known charmonium states that feed the ip' signal. 

One possible explanation for this ip' anomaly at CDF is that the ip' signal is fed by 
the decays of higher charmonium states that have not yet been discovered.0 Among 
the candidates are D-wave states, higher P-wave states, and ccg states. The main 
difficulty with this solution is to explain how these states can have such a dramatic 
effect on ip' production at the Tevatron and not have shown up in any charmonium 
experiments at lower energy. 

Another possibility is that the solution to the ip' anomaly lies in a new production 
mechanism. Braaten and Fleming have proposed that the ip'^s come primarily from 
a color-octet term in the gluon fragmentation function for -^'.0 The basis for this 
proposal is a general theory of inclusive quarkonium production that was recently 
developed by Bodwin, Braaten, and LepagcS This approach allows one to calculate 
not only perturbative corrections to any order in a^, but also relativistic corrections 
to any order in v"^. The inclusive cross section for producing a quarkonium state H 
satisfies a factorization formula of the form 

da{H + X) = Y.da{cc{n)+X) {O^), (5) 

n 

where da is the inclusive cross section for producing a cc pair separated by a distance 
less than 1/mc and in a color and angular- momentum state labelled by n. Since da 
depends only on short distances, it can be calculated using perturbative QCD. The 
nonperturbative factor {O^) in is proportional to the probability for a pointlike 
cc pair in the state n to form the bound state H. It can be defined rigorously as a 
matrix element in nonrelativistic QCD. The relative importance of the various terms 
in can be estimated using scaling rules that tell how the matrix elements scale 
with TJic and v. 

The factorization formula (|^) also applies to fragmentation functions, which have 
the general form 

D,^h{z) = E^-n(^) (^f)- (6) 

n 

For example, the process g ^ cc produces a cc pair in a color-octet ^5*1 state and 



gives rise to a term in the gluon fragmentation function of order a^: 

D,^h{z) = ^^6{l-z){O^CS,)). (7) 

This term was included in the fragmentation calculations for direct Xc production, 
because it is required for perturbative consistency. It was not included in the cal- 
culations for the S-wave states ip ^ind ip' on the grounds that the matrix element 
(Osi^Si)) is suppressed by v'^ relative to the nonperturbative factor |i?(0)p in the 
leading color-singlet term. However the term (0) may well be numerically important, 
because the leading color-singlet term is suppressed by a short-distance factor dn{z) 
that is of order a^. The CDF data on ip' production can be explained by including 
the term in the gluon fragmentation function for ip' and adjusting the matrix 
element (Cg to fit the data.0 The term (0) must also appear in the gluon 

fragmentation function for ip, and it may explain the CDF data on direct ip^s that do 
not come from ;i^c's.00 

6. Recent Developments 

Recently, there have been several significant experimental developments in char- 
monium at large pt at the Tevatron. The CDF collaboration has measured the ratio 
of the prompt cross sections for Xci and Xc2,l^ providing one more test of our under- 
standing of charmonium production mechanisms. They have also measured the cross 
sections for the bottomonium states T, T', and T"x3 This information is valuable, 
because the larger mass of the bottom quark changes the relative importance of the 
various production mechanisms. 

There have also been important developments in theory. The next-to-leading- 
order perturbative corrections to the fragmentation function (|^ have been calculated 
by Ma.0 Eventually all the relevant fragmentation functions should be calculated 
to next-to-leading order in as- Another important development is that Cho and 
Leibovich have calculated the terms in the parton cross sections for ij cc + k that 
correspond to the matrix element {0^{^Si))^ They used these results to extend 
the calculations of prompt charmonium production down to moderate pt and also to 
calculate the T production rate at the largest values of pr measured at the Tevatron. 
Ultimately, one would like to extend the theoretical predictions all the way down 
to Pt = 0. It may be necessary to include other color-octet matrix elements in the 
factorization formula (^) for the cross section at low px- It will also be necessary to 
resum the effects of soft gluons in order to get a cross section that vanishes properly 
at Pt = 0. 

In conclusion, data from the Tevatron is driving us toward a deeper understand- 
ing of the production of heavy quarkonium. The order-of-magnitude discrepancies 
between experiment and theory can be understood by introducing new production 
mechanisms that are based on a theory of inclusive quarkonium production derived 



from QCD. I believe that we are on our way to a comprehensive description of quarko- 
nium production in all high energy processes in terms of the quark mass, the QCD 
coupling constant, and a few well-defined phenomenological matrix elements. 
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